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The fluctuations exhibited by the cross-sections generated in a compound-nucleus reaction or, 
more generally, in a quantum-chaotic scattering process, when varying the excitation energy or 
another external parameter, are characterized by the width Tcorr of the cross-section correlation 
function. In 1963 Brink and Stephen [Phys. Lett. 5, 77 (1963)] proposed a method for its determi¬ 
nation by simply counting the number of maxima featured by the cross sections as function of the 
parameter under consideration. They, actually, stated that the product of the average number of 
maxima per unit energy range and Tcorr is constant in the Ercison region of strongly overlapping 
resonances. We use the analogy between the scattering formalism for compound-nucleus reactions 
and for microwave resonators to test this method experimentally with unprecedented accuracy us¬ 
ing large data sets and propose an analytical description for the regions of isolated and overlapping 
resonances. 

PACS numbers: 24.60.Ky, 05.45.Mt 


I. INTRODUCTION 

The ubiquity of quantum chaotic scattering is well rec¬ 
ognized owing to its fingerprints in diverse fields rang¬ 
ing from compound-nucleus scattering to electron 

transport through disordered mesoscopic samples m- 
Owing to the quasibound states of an open system that 
manifest themselves as resonances in the associated scat¬ 
tering processes, the cross-sections evince random fluctu¬ 
ations with respect to the energy E or some other generic 
parameter that are characterized in terms of the cross- 
section correlation functions. According to the statistical 
theory of compound-nucleus reactions [ 6 ] the coherence 
width Tcorr of the fluctuations is determined by the reso¬ 
nance widths. However, an alternative approach to effec¬ 
tively determine it from the number of maxima exhibited 
by the cross-sections as function of the excitation energy 
was introduced in Ref. [7] . It was proposed that the prod¬ 
uct of the number of maxima per unit resonance energy, 
iTjv, and the width Tcorr is constant, K^^cott = O.bbjy, 
where N is the number of statistically independent chan¬ 
nels participating in the reaction and is a constant. 
Taking into account finite range of data effects [ 8 ] yielded 
for a system with a large number of reaction channels 
boo = 2\/3l7r ^ 1.1. For two reaction channels a value 
of 62 = 0.78 was obtained [7]. The dependence of b^ 
on the experimental energy resolution of the excitation 
functions was demonstrated in Refs. laiini for a single 
channel as well as many open channels. The counting- 
of-maxima method was believed to yield more accurate 
results than the determination of Tcorr from the autocor¬ 
relation functions m- Applications to numerical exam¬ 
ples have been presented in Refs. [laiis]. 
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The main ingredient of all these findings was that one 
can infer the cross-section correlation width associated 
with a compound-nucleus reaction by simply counting 
the maxima emerging when plotting the cross-sections 
versus the excitation energy. This procedure, actually, 
was proposed at the time [7] as an add-on to the stan¬ 
dard analysis of the cross-section fluctuations [ 6 ] and was 
restricted to the Ericson region of strongly overlapping 
resonances. Recently, the counting-of-maxima method 
was extended to the regions of isolated and overlapping 
resonances pAlIlG] . There, an analytical expression for 
the dependence of the product of the number of maxima 
and the width of the associated cross-section correlation 
function as function of the energy or, more generally, of 
some parameter, on the tunneling probability was de¬ 
rived for the cases, that the latter has a Lorentzian and 
a squared Lorentzian shape, respectively. The findings 
of Refs, [m [15] actually inspired the investigations de¬ 
scribed in the present article. 

Based on the premise that the ^'-matrix formalism for 
microwave resonators m is equivalent to that originally 
derived for compound-nucleus reactions m , we exam¬ 
ined the relation of the mean density of maxima in the 
cross-section fluctuations to the correlation width Tcorr 
in chaotic scattering experiments with unprecedented 
accuracy. The experimental data, obtained from mea¬ 
surements documented in Refs. na uni in the regions 
of isolated and overlapping resonances, illustrates par¬ 
allels with studies conducted in Refs. ED ED in the 
context of nuclear physics ElESlEll. In addition, the 
effect of time-reversal invariance (T) violation, which 
has been tested in nuclear spectra [25] and for the Er¬ 
icson region in compound-nucleus reactions [26ll3T] . is 
inspected in this regard. From a theoretical perspective, 
simulations are performed using random matrix theory 
(RMT) and the Verbaarschot-Weidenmuller-Zirnbauer 
(VWZ) approach [1], accompanied by a formal appraisal 
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of chaotic scattering on networks of quantum wires [3^ . 
so-called quantum graphs, j33f[37] . In light of recent 
studies of the statistical properties of electronic transport 
in ballistic open quantum dots El, we take cognizance of 
the analogy between universal conductance fluctuations 
described by Gaussian processes in quantum dots and 
compound-nucleus fluctuations [2]. This enables us to 
analyze and extend the main results of Ref. El, which 
were originally gleaned from the S'-matrix describing elec¬ 
tron flow in quantum dots, and exhibit their applicability 
in the wider context of microwave billiards, thereby un¬ 
derlining the subtle interconnections between these seem¬ 
ingly unrelated fields. 

We applied the counting-of-maxima method to spectra, 
that were measured with a flat microwave resonator with 
the shape of a chaotic tilted-stadium billiard, and at su¬ 
perconducting conditions in a resonator with a paramet¬ 
ric shape, respectively. The aim was to verify the applica¬ 
bility of the method in the region of strongly overlapping 
resonances mil] and to verify the validity of the analyt¬ 
ical results of Refs. experimentally in the regions 

of isolated and overlapping resonances. Note, that for 
the largest achieved ratios of the resonance width and 
spacing the fluctuation properties of the ^'-matrix ele¬ 
ments are already well described by the predictions for 
the Ericson region [38]. The experimental results are 
presented in Secs. [Ill and [V| respectively. In order to cor¬ 
roborate our findings we performed RMT simulations, 
that were based on the scattering formalism introduced 
in Sec. imi Furthermore, we analyzed spectra computed 
for open quantum graphs, as outlined in Sec. IIVI 


II. FLUCTUATIONS IN AN OPEN 
MICROWAVE BILLIARD WITH AND WITHOUT 
rVIOLATION 

To realize a chaotic scattering system, a flat mi¬ 
crowave resonator with the shape of a tilted-stadium bil¬ 
liard [H uni EH was chosen, which is shown schemati¬ 
cally in the inset of Fig. The modes in the resonator 
were coupled to the exterior via two antennas denoted 
by 1 and 2, that were attached to it at the positions 
marked by crosses. For the determination of the scat¬ 
tering (S')-matrix elements S 21 describing the scattering 
process from antenna 1 to antenna 2, a vector network 
analyzer coupled microwave power into the resonator via 
antenna 1 and determined the relative phase and ampli¬ 
tude of the transmitted signal at antenna 2. The fre¬ 
quency range of 1 — 25 GHz was chosen such that only 
the component of the electric field strength perpendicular 
to the top and the bottom plate was excited. Then the 
Helmholtz equation is scalar and mathematically equiv¬ 
alent to the Schrodinger equation of the corresponding 
quantum billiard [39h4T] . Therefore, such resonators are 
called microwave billiards. The dynamics of the corre¬ 
sponding classical billiard is chaotic [42]. Thus, accord¬ 
ing to the Bohigas-Giannoni-Schmit conjecture [43], the 
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FIG. 1: Upper panel: transmission spectrum for the frequency 
interval 19 — 20 GHz, where the resonances obviously overlap. 
Lower panel: cross-section autocorrelation function C 2 i(s) for 
the same frequency interval (solid black line) and the interval 
of 13 — 14 GHz (dashed red line). Here, e gives the excitation 
frequency increment in units of the mean resonance spacing. 
The secular variation of the latter was negligible in a 1 GHz 
window. The broadening of the correlation width with in¬ 
crease in frequency is palpable. Inset: the tilted stadium 
billiard (schematic). The crosses mark the positions of the 
antennas, 1 and 2. A movable scatterer (gray) was utilized to 
obtain independent realizations. 


fluctuation properties of the eigenvalues of the associ¬ 
ated quantum billiard, i.e., of the resonance frequencies 
in the microwave billiard, are described by random matri¬ 
ces from the Gaussian orthogonal ensemble (GOE) [44] . 
We confirmed this by analyzing the resonance spectra in 
the regions of isolated and overlapping resonances m- 
An ensemble of several chaotic systems was obtained by 
introducing a scatterer (gray disk in the inset of Fig. 
into the microwave billiard and moving it to 8 different 
positions [45] . We also investigated fluctuation properties 
in chaotic scattering systems with partially broken time- 
reversal symmetry. To induce T violation we inserted a 
ferrite into the microwave billiard and magnetized it with 
an external magnetic field. 

In the upper panel of Fig. [^a measured transmission 
spectrum, i.e., the squared modulus of the ^'-matrix el¬ 
ement from antenna 1 to antenna 2, is plotted in the 
frequency range 19 — 20 GHz. The ratios of the widths 
and the spacings of the resonances depend weakly on the 
excitation frequency. At low frequencies the resonances 
were isolated, that is, their mean width was small as 
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compared to their mean spacing. However, on increasing 
it, the resonances started to overlap, as is clearly visible 
in Fig. [2 The lower panel of Fig. shows two exam¬ 
ples for the correlation function of the squared modulus 
of the S'-matrix element 5 ' 2 i, which corresponds to the 
cross section for the transition from channel 1 to channel 
2 in nuclear physics, ( 721 = | 5 ' 2 ip. It is given by [ 6 l[T3]. 

C2i{e) = {\S2,{E)\^\S2iiE + e)f) - (|52l|^)^ (1) 

where (...) denotes averaging over E. In nuclear physics 
E and E ^ s denote energies, whereas in the microwave 
experiments they correspond to the rescaled frequen¬ 
cies /, obtained by unfolding them to mean resonance 
spacing unity with the help of Weyl’s formula for the 
smooth part of the resonance density (p) of microwave 
billiards [46] which increases linearly with the frequency 
/, (p) A7rf/{2c^) with c the velocity of light and A 

the area of the billiard. In the limit of a large number 
of open channels the shape of the cross-section corre¬ 
lation function C 2 i{£) is well decribed by a Lorentzian, 

C2l{s) = C2l(0)/ + (^/Fcorr)^^- 

The averages of the resonance widths and spacings 
were observed to be approximately constant in 1 GHz 
frequency intervals. Consequently, the correlation func¬ 
tion was evaluated in a frequency range of 5 — 25 GHz 
in 1 GHz windows to ensure that the secular variation of 
the resonance parameters is negligible, yielding 10 ^ data 
points each. As illustrated in Fig. the width of the 
cross-section function increases with the frequency, that 
is, it is narrower in the frequency range 13 — 14 GHz 
(dashed red line) than for 19 — 20 GHz (solid black line). 

The number of maxima in the transmission spec¬ 

tra was also counted in 1 GHz windows. One thereby 
obtains the density of maxima, p™^^, on simply dividing 
A'max py corresponding interval length of the rescaled 
frequencies. Furthermore, each of the 1 GHz windows 
was subdivided into 10 intervals and in each of them the 
width of the cross-section autocorrelation function at half 
maximum, i.e., the value of e with C 2 i{£) = 0 . 5 C 2 i( 0 ), 
was determined in units of the mean resonance spacing. 
To further increase the data sets, transmission spectra 
were measured for 8 positions of the movable scatterer. 

Figure presents the thus obtained experimental re¬ 
sults for the product of the mean density of maxima 
and the mean correlation width (turquoise dots). The 
largest value of the latter achieved in the experiments 
was Fcorr — 1-7. We demonstrated in Ref. [38] that in 
the corresponding frequency range the ^'-matrix corre¬ 
lation functions and the distributions of the S'-matrix 
elements are well described by the predictions for the Er- 
icson region. Similarly, the product of the mean density 
of maxima and the mean correlation width approaches 
the value predicted in that limit, (p^^^)Fcorr — 0.5 600 = 
V^/tt, shown as dash-dotted line in Fig. Note, that 
(p™^^)Fcorr does uot depend on the frequency scale. 

In order to obtain an analytical expression for the de¬ 
pendence of the product (p^^^) Fcorr on the correlation 



FIG. 2: Product of the mean density of maxima and the mean 
correlation width in the transmission spectra of the tilted- 
stadium billiard shown in the inset (same as in Fig. [^ times 
the mean correlation width evaluated in 1 GHz windows in 
the T preserved case (turquoise dots). Statistical errors arise 
due to the variability of the results within the 8 realizations, 
i.e., positions of the scatterers, and are at most of the size 
of the symbols. Therefore, we omitted error bars. The red 
triangles, obtained from RMT simulations with an ensemble 
of scattering matrices of the form Eq. Q for the GOE case 
describe the experimental results well. The value of O.bboo = 
~ 0.55 (dash-dotted line) approached above Fcorr ~ 1.5 
is in good agreement with the result obtained in Ref. [7] in 
the context of nuclear physics. The dashed curve illustrates 
the analytical expression na given in Eq. § with ao = 1 and 
ai = 1/3. 


width Fcorr we used the ansatz 


/ max\ p E corr 

\Ps H corr ~ 


TT \/^lovv + 


( 2 ) 


with ao and ai as fit parameters. The comparison of our 
experimental data with the analytical expression Eq. ([^ 
revealed a very good agreement for ao = 1 and ai = 1/3. 
Hence, like in the case of the experimental data, the ana¬ 
lytical expression saturates at the value predicted for the 
Ericson region. It is shown as dashed line in Eig. [^and 
also in Eig.[^ where the product (p^^^) Fcorr is plotted for 
the case with T violation, induced by inserting a ferrite 
into the microwave billiard (black disk in EigJ^ and mag¬ 
netizing it with an external magnetic field [Si SZ] • The 
agreement of the analytical expression Eq. © and the 
experimental result for (p^^^) Fcorr again is very good for 
ao = 1 and ai = 1/3. No significant difference is notice¬ 
able between the cases of Tpreservation and T violation, 
either in the saturation value or in the general trend of 
the data points. 

In Ref. [16] the conductance fluctuations where inves¬ 
tigated in quantum dots with chaotic behaviour. There, 
an analytical expression was derived for the product of 
the mean density of maxima (p^^) and the width 7 of the 
generalized conductance correlation function which, sim¬ 
ilarly to C 2 i(e) in Eq. 0 , takes the shape of a Lorentzian 
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FIG. 3: Same as Fig. for the cases of Tviolation in the 
experiment (turquoise dots) and the GUE in the RMT simu¬ 
lations (red triangles). Inset: same as in Fig. In addition 
a magnetized ferrite (black disk) was inserted to enforce T 
violation. 


in the Ericson region, 




/ gp - i 8 r +10 

TT y 5r^ — aiT + 6 


(3) 


Here, F denotes the tunneling probability, which corre¬ 
sponds to the transmission coefficients Ta = 1 — KS'aa)!^ 
associated with antennas a = 1,2 in our microwave ex¬ 
periments, F = , where Ti T 2 . In order to verify 

Eq. we plotted in Fig. [^the ratios x(r) of the exper¬ 
imental results for (p^^^)rcorr for the cases of preserved 
T (turquoise dots) and of T violation (red triangles) and 
the prediction Eq. x(r) = (pf*'’')rcorr/((/3?°)7)- For 
F > 0.4 the ratio takes the value x(r) 1, i.e., there the 
agreement between the experimental result and the pre¬ 
diction Eq. ^ is very good. For small values of F devi¬ 
ations are clearly visible and, actually, expected because 
for F —> 0 the product {p^^^)Tcorr vanishes, see Figs. 
and[^ whereas the square-root function in Eq. <© takes 
a non-zero value for F ^ 0 El- 


Ill. FLUCTUATIONS IN A RANDOM 
5^-MATRIX MODEL FOR CHAOTIC 
SCATTERING PROCESSES 

To corroborate our experimental findings we performed 
RMT simulations similar to those presented in Refs. m 
[20] . For this we used the ^'-matrix formalism developed 
in m in the context of compound-nucleus reaction the¬ 
ory, which, actually, also has been used in [HUS] for the 
RMT simulations. The associated unitary S matrix has 
the general form 



FIG. 4: (Golor on line) Ratio x(r) = (p™)rcorr/((p?^)7) 
of the experimental result (p™^^)Fcorr for the cases of 
T preservation and T violation (turquoise dots and green 
squares, respectively) and the analytical result for a quan¬ 
tum dot, see Eq. Eurthermore, the ratio x(r) = 

{px^^)Xc/{{Px^)Xc) of the RMT simulation for the param¬ 
eter dependent system {px^^)Xc and the analytical result 
Eq. ( p^ for a quantum dot is plotted (red triangles). Here, F 
denotes the average of the transmission coefficients associated 
with the two antennas 1 and 2, F = (Ti +T2)/2, in the case of 
the microwave billiard and the tunneling probability in that 
of a quantum dot. 


in the cavity [48]. They were chosen as real Gaussian 
distributed random numbers with zero mean. As stated 
above, the fluctuation properties of the resonance fre¬ 
quencies of the resonator are well described by random 
GOE matrices in the T-preserving case. Thus we in¬ 
serted in Eq. Q for the matrix which actually stands 
for the Hamiltonian of the closed resonator, a real and 
symmetric random matrix of dimension N = 200 from 
the GOE [ 44 ]. In order to study a chaotic scattering 
system with induced partial T violation we proceeded as 
in Refs. da na and inserted for H in Eq. @ a com¬ 
plex random matrix of the form H = + 

with a member from the GOE, a real, antisym- 
metric matrix with Gaussian-distributed entries and ^ 
the T -breaking parameter. The number of open chan¬ 
nels M, that is, the dimension of S{f) was chosen equal 
to 32, i.e., we introduced 30 fictitious ones. Further¬ 
more, since in the experiments the average S matrix was 
diagonal, we required that this is also the case for the 
matrix WW^ dJiiHi. Its entries are solely determined 
by the transmission coefficients Tc, with c = 1, • • • , M. 
For the RMT simulations we used the same values of the 
transmission coefficients as in the experiments within the 
1 GHz frequency intervals, see Ref. [20] . 


S{f) = l-2mW{fl-H+ iTrW'^W) . (4) 

The matrix elements of W describe the coupling of the 
antennas to the resonator modes 113 and comprise the 
fictitious channels that account for the Ohmic absorption 


The RMT results plotted as red triangles in Figs. 
and [ 3 ] were obtained with an ensemble of 1000 random 
scattering matrices of the form Eq. 0. The data pre¬ 
sented in Fig. were generated by choosing ^ as in the 
measurements uniEni, where it varied with the excita- 
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tion frequency. The RMT results for ^ = 0 and ^ 7 ^ 0 
are indistinguishable. Thus, the value of the product 
^^max)r^orr docs uot dcpcud ou whether Tis preserved 
or not, in agreement with the experimental results. It 
saturates at 0.5 600 = i.e., at the value obtained 

in the Ericson region in the context of nuclear physics 
for systems with a large number of reaction channels or 
open channels [ 8 ]. Furthermore, the RMT results are 
again well described by the formula Eq. with ao = 1 
and ai = 1/3. 


IV. FLUCTUATIONS IN OPEN QUANTUM 
GRAPHS 


As a second test system we analyzed the S matrix of 
open quantum graphs, networks of one-dimensional wires 
(bonds) joined at vertices [32] and opened by attaching 
leads that extend to infinity at some of the vertices. An 
example is shown in the inset of Fig. where a so-called 
tetrahedral graph is presented, which consists of 6 bonds 
B, 4 vertices V and has a lead L attached to the lower 
left vertix and to the center one, respectively. Quantum 



FIG. 5: Product of the mean density of maxima and the 
mean correlation width for a chaotic tetrahedral graph shown 
schematically in the inset, which consists of 6 bonds B, 4 ver¬ 
tices V and has two leads L attached, corresponding to two 
open channels, obtained by varying the coupling w between 
the leads and the vertices. Shown are the results for the case of 
T preservation (turquoise dots) and for complete T violation 
(red squares), where for all bonds Aij = 0.5, i < j was used. 
The data are well described by Eq. Q with do = 1, di = 1/3 
and d 2 = 1/6 (dashed line). For large Fcorr the product ap¬ 
proaches the value 0.5^2 = 0.39 (dash-dotted line) predicted 
for nuclear reactions with 2 independent channels [7]. Note 
that the graph shown in the inset is merely representational 
since the lengths do not reflect those used in the actual sim¬ 
ulations. 

graphs can be handled mathematically and at the time 
serve as a simple model for chaotic scattering. They, in 
fact, have the particular property that the semiclassical 
approximation is exact [49] and, moreover, the S'-matrix 
correlation functions could be calculated explicitely and 
shown to agree with the RMT result [T] [35j [37l [50] . The 
mathematical model for the S matrix of open quantum 


graphs has been formulated in Refs. [331EI]. The com¬ 
ponent 2 pij of the graph wavefunction 7 / for the mode 
propagation in the bond from vertex i to vertex j satis¬ 
fies the one-dimensional Schrodinger equation 

/ ^ \ 2 

f “ (5) 


where Aij = —Aji is a magnetic vector potential with a 
nonvanishing real part introduced on the bonds to induce 
T violation. We set it to zero for the study of chaotic 
scattering in T -preserved systems, and chose it purely 
real with \Aij\ being the same on all bonds to induce 
Tviolation, otherwise. Furthermore, we required that 
the wave functions are continuous and satisfy Neumann 
boundary conditions at the vertices in order to ensure 
current conservation. On attaching a lead to M of the in 
total V vertices of a graph, the S matrix can be expressed 
as 

S{f) = t-2TTiW{h{f) + iTTW^Wy'^W^, ( 6 ) 

with W being the M x V leads-vertices coupling matrix 
introduced in Ref. [33] for a tunable coupling [52], i.e., its 
matrix elements depend on a parameter re, 0 < Wi^j = 
dij^jw/y/iT < 1 , where iL denotes the indices of the ver¬ 
tices with leads. The term h{f) is the Hamiltonian of the 
closed graph with the matrix elements 


%(/) 


{ ^m^i Cim COt ( Lim^ ? ^ — j 

Cij sin“^ (^-^ij) > ^^3 


Here, C is the connectivity matrix with entries Cij = 1 
if the vertices i and j are connected, and zero otherwise. 
Furthermore, L is the length matrix of the graph, where 
Lij gives the length of the bond joining the vertices i and 
j. Note the similarity of the form of the S matrix of an 
open graph, Eq. (§, with that of a compound-nucleus 
scattering process given in Eq. the only difference 
being their dependence on the frequency through the in¬ 
verse operator. In the former case it coincides with that 
of the graph Hamiltonian h(/), whereas in the latter case 
it depends linearly on the frequency, while the Hamilto¬ 
nian of the closed system is modelled by a frequency- 
independent random matrix. 

We considered in our numerical simulations tetrahedral 
graphs consisting of 4 vertices and graphs that comprised 
60 vertices, where in both cases the connectivity was 3 
(see insets of Figs. [^and[^. The lengths of the bonds 
were confirmed to be rationally independent by choosing 
values which were the square roots of prime numbers. In 
addition, the chaoticity of the scattering dynamics was 
independently verified by the fluctuation properties of 
the eigenvalues of the closed graphs. The ^'-matrix ele¬ 
ments were computed over a large frequency range and 
the frequencies were again rescaled to mean resonance 
spacing one using WeyFs formula for the mean resonance 
density (p) in quantum graphs, which is frequency inde¬ 
pendent, (p) >C/(27r), with C the total length of the 
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graph. Consequently the mean resonance spacing is the 
same for the whole spectrum, in stark contrast to the sit¬ 
uation in quantum (microwave) billiards. The same holds 
for the correlation width. Accordingly an ensemble of 
statistically indendent data sets was obtained by simply 
dividing the spectrum into frequency windows containing 
a sufficiently large number of resonances, determining in 
each of them the number of maxima and the mean width 
of the cross-section correlation function and then calcu¬ 
lating their arithmetic averages. 

In the case of the tetrahedral graph shown schemati¬ 
cally in the inset of Fig. we attached to two vertices a 
lead and varied the size of Fcorr by altering the parameter 
w of the entries of the W matrix, which physically corre¬ 
sponds to changing the strength of the coupling between 
the leads and the vertices or that of the graph to its en¬ 
vironment. In Fig. the quantity (p^^^)rcorr is plotted 
against Fcorr for the T-preserving (turquoise diamonds) 
and the T -violated (red squares) cases. Especially for 
the T -invariant graph its saturation at 0.5 62 ~ 0.39 
(dash-dotted line) is clearly visible. This value was pre¬ 
dicted for 2 reaction channels in nuclear physics [7] and 
for 2 — 3 open channels in quantum dots M- Again, 
the product is well described by an expression similar to 
Eq. 


{pTn^c 


0.39do 


y/r^orr + ^^iTcorr + 


( 8 ) 


The fit of Eq. <§ to the numerical data yielded do = 
1, di = —1/3, ^2 = 1/6 (dashed line). Thus, the analyt¬ 
ical expression saturates at the same value, 0 . 662 , as the 
experimental data. 

In order to attain higher values of Fcorr, we used the 
fact that it increases with the number of vertices V [34] . 
Accordingly, we considerered a graph with V = 60 which 
also accorded the flexibility of increasing the number of 
open channels. The product Fcorr was calculated 

for different values of the correlation width by attaching 
leads at random to a varying number of 2—15 vertices and 
also altering the value of w in the elements of the W ma¬ 
trix. Figure which compares the resulting values with 
the RMT simulations already presented in Fig. reiter¬ 
ates the agreement with Eq. (§ with ao = 1, ai = 1/3 
shown as dashed line and, furthermore, with the results 
obtained in nuclear physics a few decades earlier (dash- 
dotted line) for Fcorr ^ 1-75, and highlights the anal¬ 
ogy between open channels in generic chaotic scattering 
systems and the reaction channels in compound-nucleus 
scattering processes. The results presented in this sec¬ 
tion were obtained by analyzing numerically obtained S- 
matrix elements. Currently, we build up an experimental 
setup for high-precision measurements of the ^'-matrix 
elements of quantum graphs using networks of supercon¬ 
ducting waveguides. 

In the previous sections 021 and[lV|we investigated 
the fluctuations exhibited by the cross sections when 
varying the excitation frequency. In the following sec¬ 



FIG. 6 : Product of the mean density of maxima and the 
mean correlation width obtained from simulations performed 
for a chaotic graph with y = 60, shown schematically in the 
inset, and a varying number of open channels (2 — 15) and 
couplings w between the leads and the vertices are illustrated 
as turquoise dots. The red triangles are the same as those 
exhibited in Fig.[^ A saturation at 0.5 600 = 0.55 is observed, 
in agreement with the value cited in Ref. [ 8 ]. Furthermore, 
the data are well described by Eq. § with ao = 1, ui = 1/3. 
Note that the equality of the side-lengths of the graph shown 
in the inset is solely for representational convenience. 


tion we study the counting-of-maxima method in a sys¬ 
tem that depends on a geometric parameter. 


V. FLUCTUATIONS IN AN OPEN 
PARAMETRIC MICROWAVE BILLIARD 

The analysis of the fluctuation properties of the eigen¬ 
values of closed, parameter dependent systems provided 
new insight into the universal behavior of chaotic sys¬ 
tems [21 [53ll56] . In order to study cross-section fluctua¬ 
tions of an open system with respect to a parameter, as 
opposed to the frequency, we analyzed the spectra mea¬ 
sured with a microwave billiard having the shape of a 
desymmetrized straight-cut circle presented in Ref. m 
and shown schematically in Fig.[^ the classical dynamics 
of which is chaotic ISIIISH]. The varying parameter, a, 
concured with the angle of rotation of a dielectric wedge 
of Teflon® inside the resonator and assumed 37 equally 
spaced values in steps of 2.5°. In order to ensure a suf¬ 
ficiently high spectral resolution and to eliminate dissi¬ 
pative processes, the experiment was performed at 4.2 K 
with a superconducting microwave billiard manufactured 
from lead-plated copper m- Thus, complete sequences 
of 440 eigenfrequencies could be identified for each value 
of a. The fluctuating positions of the resonance frequen¬ 
cies trace out irregular oscillatory curves that avoid in¬ 
tersections as the parameter a is varied (see Fig. 2 of 
Ref. [57]). 

The eigenfrequencies were unfolded to mean spacing 
unity for a fixed value of a with the help of Weyl’s for¬ 
mula, yielding 440 rescaled eigenfrequencies ei{a), i = 
1 , 2 , • • • , 440 for each of the 37 values of a. Also the pa- 
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FIG. 7: Product of the mean density of maxima and the 
mean correlation width obtained for the parameter-dependent 
cut-circle system (turquoise dots) at 440 different frequencies. 
The statistical errors are no larger than the size of the sym¬ 
bols. The dashed curve illustrates the analytical expression 
Eq. 0 for the variation of an external parameter which de¬ 
scribes the data well, with co = 1, ci = —1/2 and C 2 = 1/3. 
Upper inset: sketch of the billiard along with the rotatable 
wedge-shaped piece of Teflon®. The angle of rotation of the 
wedge dehnes the parameter but its initial orientation is ar¬ 
bitrary. Lower inset: zoom into the experimental data range 
Xc < 0.08. 


rameter a was rescaled by using the procedure described 
in Ref. [59] , 



We considered the cross-section autocorrelation function 
for parametric variations, 

C^^{5X) = {\S2i{X)\^\S2i{X + 8X)\^) - (|52iH2, (10) 

where (...) denotes the average over the parameter and 
intermediate frequencies between the respective eigen¬ 
value and the neighboring ones. This correlation func¬ 
tion has been examined for quantum dots and has 
a square Lorentzian shape 0 na na, c^i{5x) = 

C^(0)/ ^1 + {5X/Xc)^^ . The correlation width Xc of 

C 2 i(SX) was either obtained from a fit of this formula to 
it or by determining the value of 6X for which C^i(SX) = 
C^(0)/4. Both procedures yielded the same values of 
Xc- Maxima were counted in the spectrum \S 2 i{X)\‘^ for 
each frequency, using the fact that the eigenfrequencies 
correspond to the frequencies at the maxima in the res¬ 
onance spectra, |5'2i(/)P- Accordingly, we determined 
the number of maxima in a given X range by simply 
counting the eigenfrequencies, which actually had been 
determined for the studies published in Ref. m, and 
the density of maxima, was inferred on dividing 

the number thereof by the range of the rescaled parame¬ 
ter. Figure]^ shows the quantity {px^^)Xc as a function 
of the correlation width Xc- It is expected to assume 
a constant value of 3/(7r\/2) 0.68 (dash-dotted line) 


for Xc ^ 1. The absence of such saturation may be at¬ 
tributed to the attainment of only very small values of 
Xc, since the measurements were performed at supercon¬ 
ducting conditions. Consequently, there is no dissipation 
in the walls, implying that the effective number of open 
channels is 2. This is in contrast to the experiments with 
the tilted-stadium microwave billiard at room tempera¬ 
ture, where fictitious open channels had to be introduced 
to account for the absorption. For very small correlation 
widths Xc an analytical expression, which is similar to 
that given in Eq. (§ was derived in Ref. [15], 

To test the results for the parametric system in more 
detail, we in addition performed RMT simulations using 
the S'-matrix formalism Eq. We replaced by a 

Hamiltonian H{p) = Hi cos p + H 2 sin/i, where iLi, H 2 
are random matrices drawn from the GOE and p is var¬ 
ied m- The number of open channels was chosen equal 
to M = 32. Still the results of these simulations, plotted 
as red triangles in Eig. describe the experimental data 
very well and saturate for Xc ^ 0.7 at {px^^)Xc — 0.68, 
i.e., at the value predicted for the Ericson region. Like 
in the experiments with the tilted stadium in Sec. HII we 
derived an analytical expression for the dependence of 
the product {px^^)Xc on the correlation width Xc in 
the regions of isolated and overlapping resonances using 
a similar ansatz 

max\ Y' _ ^ 

^ ^Xl + ciXc + C2 

with now 3 fit parameters cq, ci, C 2 - The good agree¬ 
ment with Eq. (ITh (dashed line), obtained by setting 
Co = 1, Cl = — 1/27 and C 3 = 1/3 is clearly visible. Eor 
large values of Xc, the RMT results for {p^'^)Xc and 
the analytical expression attain the predicted asymptotic 
value, l\p^'^)Xc experimental data are also 

well described by this expression. This is further demon¬ 
strated in the inset which shows a zoom into the experi¬ 
mental data range. 

In Refs, [m |T6] an analytical expression was derived 
for the product of the mean density of maxima (p^) and 
the width Xc of the correlation function for the variation 
of an external perpendicular magnetic field in quantum 
dots, which takes the shape of a square Lorentzian in the 
Ericson region. 


{pT)Xc 


V3 /7r2 - ior + 6 
2r2 - 3r + 2 ’ 


( 12 ) 


with F the tunneling probability, which corresponds to 
the average of the transmission coefficients Ti and T 2 as¬ 
sociated with the antennas in the microwave experiment. 
Figure 1^ shows the ratio x(F) = {p^^^)Xc/{{p%^)Xc) 
of the experimental result and the analytical prediction 
Eq. ( [T^ for the product of the mean density of maxima 
and the mean correlation width. It illustrates that, like 
in the parameter-independent cases considered in Sec. El 
the former are well described by the latter for F > 0.4. 
While for quantum dots, the product of the mean density 

















of maxima and the width of the corresponding correlation 
function attains a non-vanishing value for F ^0, this is 
not the case for an open microwave or quantum billiard. 
Therefore, the deviation of x(r) from unity, observed for 
small values of F, was expected. 

VI. CONCLUSIONS 

Drawing inspiration from ideas developed in nuclear 
physics, in this article the relation of the mean den¬ 
sity of maxima in cross-section fluctuations of a 

chaotic compound-nucleus scattering process to the cor¬ 
relation width Fcorr using the analogy of the associated 
^'-matrix formalism to that for the description of the 
resonance spectra of microwave billiards has been inves¬ 
tigated. Historically, the value of = 2iFArFcorr has 
been extensively debated [ 31 ] in the nuclear physics lit¬ 
erature. Therefore it is reassuring that our experiments 
were able to reconfirm the prognostication of Refs. laiHi 
regarding the saturation value for fairly large values of 
Fcorr that were predicted nearly half a century ago. Not 
only do we establish the association thereof with conduc¬ 
tance fluctuations in quantum dots but our numerical 
and RMT results also suggest an appropriate ansatz to 
account for the behavior in the region of low Fcorr- In¬ 
deed, on the basis of the fit of a general ansatz for the 
product of the mean density of maxima and the mean 
correlation width to the experimental and RMT results 
a good analytical description was found, which inter¬ 
polates between zero for a vanishing correlation width 
and the predicted saturation value in the Ericson region, 
for the cases of a large number of open channels and 


of just 2 open channels and for a parameter dependent 
system, see Eqs. (§, (§ and & respectively. Interest¬ 
ingly, all three expressions have the same simple struc¬ 
ture and also the coefficients are similar. Still, we were 
not able to provide an analytical derivation. We also 
compared our experimental and numerical results with 
predictions for quantum dots m associated with a cor¬ 
relation function with the shape of a Lorentzian and a 
squared Lorentzian, respectively, and found a good agree¬ 
ment for sufficiently large tunneling probabilities, i.e., 
transmission coefficients F ^0.4. For small F -> 0, how¬ 
ever, the product of the mean density of maxima and the 
mean correlation width vanishes in our experiments and 
RMT simulations and takes a non-zero value in quantum 
dots. The extension of the results to quantum graphs 
further broadens their scope of applicability to studies 
in condensed matter and solid-state physics for which 
quantum graphs have gained widespread acceptance as 
model systems [60]. Hence, we believe that the results 
presented herein hold great potential at the interface of 
multiple fields for future investigations in this direction, 
while simultaneously bringing long-awaited closure to a 
lingering problem in nuclear physics. 
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